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In this paper we review a recently developed approximate method for investigation of 
dynamics of compressible ellipsoidal figures. Collapse and subsequent behaviour are de¬ 
scribed by a system of ordinary differential equations for time evolution of semi-axes of 
a uniformly rotating, three-axis, uniform-density ellipsoid. First, we apply this approach 
to investigate dynamic stability of non-spherical bodies. We solve the equations that 
describe, in a simplified way, the Newtonian dynamics of a self-gravitating non-rotating 
spheroidal body. We find that, after loss of stability, a contraction to a singularity occurs 
only in a pure spherical collapse, and deviations from spherical symmetry prevent the 
contraction to the singularity through a stabilizing action of nonlinear non-spherical oscil¬ 
lations. The development of instability leads to the formation of a regularly or chaotically 
oscillating body, in which dynamical motion prevents the formation of the singularity. 
We find regions of chaotic and regular pulsations by constructing a Poincar diagram. 
A real collapse occurs after damping of the oscillations because of energy losses, shock 
wave formation or viscosity. We use our approach to investigate approximately the first 
stages of collapse during the large scale structure formation. The theory of this process 
started from ideas of Ya. B. Zeldovich, concerning the formation of strongly non-spherical 
structures during nonlinear stages of the development of gravitational instability, known 
as ’Zeldovich’s pancakes’. In this paper the collapse of non-collisional dark matter and 
the formation of pancake structures are investigated approximately. Violent relaxation, 
mass and angular momentum losses are taken into account phenomenologically. We es¬ 
timate an emission of very long gravitational waves during the collapse, and discuss the 
possibility of gravitational lensing and polarization of the cosmic microwave background 
by these waves. 


1. Introduction 

Observed Universe is non-homogeneous at different scales. Perturbations of density lead 
to gravitational collapse of matter and formation of structures of different sizes: stars, 
galaxies, large scale structure. Many of astropliysical objects are essentially non-spherical, 
and their formation and evolution can be approximately described as a behaviour of self- 
gravitating ellipsoidal bodies in the frame of Newtonian gravity. 

Theory of ellipsoidal figures was developed in classical works of Newton, Maclaurin, 
Jacobi, Poincare, Chandrasechar. T he systemati c analysis of ellipsoidal figures of equilib¬ 
rium is presented in the book of 1 Chandrasekhar! ( 19691) . Different types of incompressible 


f Email address for correspondence: gkogan@iki.rssi.ru, tsupko@iki.rssi.ru 







2 G. S. Bisnovatyi-Kogan and O. Yu. Tsupko 

ellipsoids are investigated there on a basis of virial equations, as well as their dynamical 
and secular stability. 

The shape of ellipsoid is determined by sizes of its three semi-axes. Non-rotating figures 
always take the spherical form, all semi-axes are equal. In the presence of rotation the 
spheroid (axially-symmetrical ellipsoid with two equal semi-axes) is formed (Maclaurin 
spheroid). If the angular momentum of rotation is large enough, development of different 
kinds of instabilities lead to formation of a three-axis ellipsoid, with rotation around the 
small axis (Jacobi ellipsoid). Incompressible spheroid can be secularly or dynamically 
unstable for transition to three-axis ellipsoid. The points of the onset of the bar-mode 
dynamical, and secular instabilities of the Ma claurin sp hero id, for its transition into the 
Jacobi ellipsoid, are found in the book [Chandrasekhar! (1969:. 

Approach based on ellipsoidal description allows to investigate the time evolution of 
a body considering behaviour of its semi-axes. First of all, its initial collapse can be de¬ 
scribed, which is a central problem in formation of astrop hysical objects. In the ca se of 
the collapse compressible ellipsoidal objects are considered. Lvnden-Bell ( 1964 1965h ap¬ 
plied the Chandrasechar’s virial tensor method for the rotating, self-gravitating spheroid 
of pressure-free gas, and had shown the growth of non-axisymmetric perturbations dur¬ 
ing the collapse. It was also discussed that the slowly shrinking Maclaurin spheroid will 
enter the Jacobi seri es if it shrinks slowly enough for the dissipative mechanisms to be 
operative. [Ze ldovi chl (119641) considered the dynamics of rotating dust ellipsoids with lin¬ 
ear dependence of velocities on coordinates, when the problem is solved e xact ly. T here 
are nu meri cal investigations of collap sing pressureles s spheroids in paper s of Lvnden-Bell 
( 19641 ) and Lin. Mestel fc Shu ( 19651 ) . In the paper of Rosensteel fc Tranl ( 1991 ) the virial 
equations for rotating Riemann ellipsoids of incompressible fluid are demonstrated to 
form a Hamiltonian dynamic al system. There is a d e tailed descript ion of ellipsoid model 
of rotating stars in papers of lLai. Rasio fc Shapirol ( 1993 . 1994al lbl). Using a variational 
principle, they had derived and investigated equations of the evolution of a compressible 
Riemann-S ellipsoid, incorporating viscous dissipation and gravitational radiation. So¬ 
lutions of these approximate equations permitted to obtain equilibrium models, and to 
investigate their dynamical and secular stability. The point o f ins tability o f comp ressible 


investigate their dynamical and secular stabi lity, The point p t ins tability o l comp ressibl 
Maclaurin spheroids is found in Lai. Rasio fc Shapiro dl993) and Shapirol (_2004). Them 
is a wide ranging review of this topic in the paper of Lvnden-Belll ( 19961) in memory of 
Chandrasekhar. 

In the present paper we review the developed approximate model for consideration 
of dynamics of collapsing ellipsoidal figures. Present review is based on the papers of 
Bisnovatvi-Kogan ( 2004 ). Bisnovatvi-Kogan fe Tsupkol ( 2005 ). Bisnovatvi-Kogan fe Tsupko 
( 2008a)). We derive and solve equations of the dynamical behavior in a simple model of 
a compressible uniformly rotating three-axis ellipsoid. The ordinary differential equa¬ 
tions for ellipsoid semi-axes evolution in time are obtained with help of variation of the 
Lagrange function of ellipsoid. In our model a motion along three axes takes place in 
the gravitational field of a uniform density ellipsoid, with account of rotation, repre¬ 
sented by centrifugial forces, and the isotropic pressure, represented in an approximate 
non-differential way. 

Our model can be used for investigation of star dynamics, and also for the first stages 
of formation of large scale structure consisted from dark matter, with appropriate ap¬ 
proximations. Starting from section [2] we apply our model for investigation of dynamic 
stability of non-spherical stars. We show, that deviation from the spherical symmetry 
in a non-rotating star with zero angular momentum leads to formation of a regularly or 
chaotically oscillating body, in which dynamical motion prevents the formation of the 
singularity. The non-spherical star without dissipative processes never will reach a singu- 
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larity in the Newtonian dynamics. Therefore collapse to a singularity is connected with 
a secular type of instability, even without rotation. Starting from section 0 we consider 
approximate model for large scale structure formation. We estimate emission of a very 
long gravitational waves during a collapse of the dark matter, and discuss a possibility 
of gravitational t ensi ng, a nd p olarisation of CMB by these waves. 

Recently lYandervoortl ( 201lj ) models the galaxy as a heterogeneous ellipsoid with an 
arbitrary stratification of the density distribution, and represent the mean motions of 
the stars in terms of a velocity field that sustains that density distribution consistently 
with the equation of continuity. The virial equations are reduced to a closed system of 
equations of motion that govern the semi-axes of the characteristic ellipsoid. The semi¬ 
axes of an ellipsoid characterizing the size and shape of the galaxy are functions of the 
ti me, w hereas the stratification of the density does not change. In the second paper 
( VandervoortJ 12014 1 author describes an application of the model to an investigation 
of chaotic behaviour in axisymmetric pulsations of non-rotating, spherically symmetric 
config urati ons a nd of rotat in g, axisym m etric c o nfig u rations . See a lso papers of R odrig ues 
( 2014) . Sharif fe Zaeem U1 Haa Bhatti ( 2013 1. Sharif &; Zaeem U1 Haa Bhattil (|2014f) . 


2. Dynamic stability of non-spherical stars 

Dynamic stability of spherical stars is determined by an average adiabatic power 

<91ogP 


7 = 


dlogp 


For a density distribution 


P = PoP 


(£)■ 


( 2 . 1 ) 

( 2 . 2 ) 


the star in the newtonian grav ity is stable against dynamical collapse when (1 Zeldovich &; Novikov 
1967 : Bisnovatvi-Kogan 19891) 



P 


dm. 




> 0 . 


(2.3) 


Here po is a central density, M is a stellar mass, m is the mass inside a Lagrangian 
radius r, so that m = 4n / Q r pr 2 dr, M = m(R ), R is a stellar radius, ip (^) is the Emden 
function corresponding to 7 = 4/3. 

Let us consider equation of state P = I\ p 1 . Spherical star with such equation of state 
is stable against collapse if 7 > 4/3 and unstable at 7 < 4/3, when the instability leads 
to a collapse to singularity. At 7 = 4/3 there is only one equilibrium state for a given 
mass, at a unique value of the constant A'o. This state is neutral to the stability, and 
may suffer contraction to singularity at K < A'o, or expansion to infinity at AT > Kq. 

Collapse of a spherical, initially unstable star may be stopped only by a stiffening of 
the equation of state, like neutron star formation at late stages of evolution, or formation 
of fully ionized stellar core with 7 = | at the collapse of clouds during star formation. 
Without such stiffening a spherical star in the newtonian theory would collapse into a 
point with an infinite density (singularity). 

In the presence of a rotation a star is becoming more dynamically stable against 
collapse. Due to the more rapid increase of a centrifugal force during contraction, in 
comparison with the newtonian gravitational force, collapse of a rotating star will be 
always stopped at finite density by centrifugal forces. 

Here we show, that deviation from the spherical symmetry in a non-rotating star with 
zero angular momentum leads to a similar stabilization, and non-spherical star without 









































4 


G. S. Bisnovatyi-Kogan and O. Yu. Tsupko 


dissipative processes never will reach a singularity. Therefore collapse to a singularity is 
connected with a secular type of instability, even without rotation. 

We calculate a dynamical behavior of a non-spherical, non-rotating star after its loss 
of a linear stability, and investigate nonlinear stages of contraction. We use approximate 
system of dynamic equations, d escribing 3 degrees of freedom of a un iform self-gravitating 
compressible ellipsoidal body ( Bisnovatvi-Koean fc Tsuokol l2008al) . We obtain that the 
development of instability leads to the formation of a regularly or chaotically oscillat¬ 
ing body, in which dynamical motion prevents the formation of the singularity. We find 
regions of chaotic and regular pulsations by constructing a Poincare diagram for dif¬ 
ferent values of the initial eccentricity and initial entropy. For simplicity we restrict 
ourself by calculating only spheroid al figures with 7 = |, see also results for 7 = | in 
iBisnovatvi-Kogan fc Tsupkol ( 2008all . 


3. Equations of motion for semi-axes of rotating 3-axis ellipsoid 

Here we derive equations of motion for general case of rotating 3-axis ellipsoid and 
consider cases 7 = 5/3 and 7 = 4/3. 

Let us consider a compressible 3-axis ellipsoid with semi-axes a ^ b ^ c 



(3.1) 


We assume that the ellipsoid rotates uniformly with an angular velocity SI around the axis 
z. Let us approximate the density of the matter p in the ellipsoid as a uniform. During 
collapse and posterior time evolution the sizes of semi-axes a, b , c, angular velocity 
S7 and density p are changed with time, but ellipsoid keeps ellipsoidal shape with the 
homogeneous distribution of density and uniform rotation. 

The constant mass m and total angular momentum M of the uniform ellipsoid are 
connected with uniform density, angular velocity and semi-axes as (V is the volume of 
the ellipsoid) 


m = pV 


4 "7r C 77~i 

— pabc , M = pQ ( x 2 + y 2 ) dV = — H(a 2 + b 2 ). 
v 


(3.2) 


We assume a linear dependence of the velocities on the coordinates in the rotating frame 


by 


(3.3) 


If an homogeneous ellipsoidal body of inviscid fluid initially has velocities that are linear 
functions of position, then it will remain ellipsoidal and the internal velocities will re¬ 
main linear functions of p ositio n for all time provided that the density remains spatially 
uniform (see, for example, Lvnden-Belll ( 19961 V). 

In absence of any dissipation this ellipsoid is a conservative system. To derive the equa¬ 
tions of motion (equations describing time evolutions of semi-axes) let write a Lagrange 
function of the ellipsoid which consists of kinetic energy Ukin and potential energy U po t 
of ellipsoid: 


U — Ukin Up 0 t > Upot — Ug -f- Eth T U ro t ■ 


(3.4) 


The kinetic energy is written as (with using (|3.3[0 


Ukin — 2 P 


(vl +vl + ' 


l)dV = ^(c 


+ b 2 +c 2 ). 


V 


(3.5) 
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The gravitational energy U g of the uniform ellipsoid is defined as ( Landau fc Lifshit3 
19931: 


Ug = - 


3Gm 2 


du 


10 J yj (a 2 + u){b 2 + u){c 2 + u) 


(3.6) 


In case of spheroids (a = b c) this in tegral can be taken analy ti cally, in case of ellipsoid 
it is ex pressed in elliptical integrals (see Bisnovatvi-Kogan ( 2004 1. Bisnovatvi-Kogan fc Tsupkol 

(f2005h V 

The rotational energy is 

Urot = \p J v rot dv = 7 }Ptt 2 J{x 2 + y 2 ) dx dy dz = ^ fl 2 (a 2 + b 2 ) . (3.7) 

v v 

The angular velocity is changed with time, and it is more convenient to express (157b via 
M which is constant in absence of dissipation. Taking into account the second expression 
in (1371) . we obtain the relation 

* M 2 


U r ot — 


(3.8) 


2 m(a 2 + b 2 ) 

Different form of thermal energy E t h are used depending on problem under consideration. 


3.1. Rotating ellipsoid, case 7 = 5/3 

First, we use the non-relativistic equation of state P = Kp 1 with polytropic index 7 = 
5/3. It corresponds to the case of a matter consisted of non-collisional non-relativistic 
dark matter particles. 

The total thermal energy, with adiabatic index n = 1.5 (7 = 5/3), of non-relativistic 
dark matter particles in the ellipsoid is Eth ~ V~ 2 / 3 ~ ( abc)~ 2 / 3 , and the relation 
between pressure P and thermal energy Eth is Eth = § PV. Let us write E t h using initial 

Values Eth,ini ®ini bini Gn a.S 


E t h 


Eth,in 


m'-'in'-'in 


) 2/3 


{abc) 2 / 3 


£ 

0 abc ) 2 / 3 ' 


(3.9) 


We introduce here e, which we will call as the entropy function. The entropy function e 
is constant in the conservative case, but it is variable in presence of a dissipation. 

By variation of the Lagrange function we obtain the Lagrange equations of motion: 


3Gm 


a = — 


2 


b 


3Gm 

2 



du 

10 1 £ 

25M 2 

a 

(a 2 + u) A 

3 m a {abc) 2 / 3 

2 

m z 

{a 2 + b 2 ) 2 

du 

10 1 £ 

25 M 2 

b 

(b 2 +u)A 

3to b {abc) 2 / 3 

2 

m z 

{a 2 + b 2 ) 2 


(3.10) 

(3.11) 


OO 


3 Gm f du 
~^~ C J (c 2 + u)A 


10 1 £ 

3 m c ( abc ) 2 / 3 ’ 


(a 2 + u)(b 2 + u)(c 2 + u). (3.12) 


0 

At the right-hand side of these equations there are forces per unit ellipsoid mass. The 
first term is the gravitational force directed to the center of ellipsoid, the second term 








































6 


G. S. Bisnovatyi-Kogan and O. Yu. Tsupko 


is the force of thermal pressure directed outward, and the third term is the centrifugal 
force directed outward. 

It is easy to check that equilibrium solution of these equations are the Maclaurin sphe¬ 
roid and the Jacobi ellipsoid ( Chandrasekhar 19691 ). 


3.2. Non-rotating ellipsoid, 7 = 4/3 

For investigation of dynamic stability of non-spherical bodies we consider the non¬ 
rotating compressible homogeneous ellipsoid with the equation of state P = Kp 7 , 7 = 
4/3. The case 7 = 5/3 considered above is not interesting for investigations of dynamic 
stability, because isentropic spherical star with 7 = 5/3 always stops contraction, and 
never suffers collapse to singularity. 

A spherical star with 7 = 4/3 collapses to singularity at small enough K , and we 
show here, how deviations from a spherical form prevent formation of any singularity. 
For 7 = 4/3, the thermal energy of the ellipsoid is E t h ~ V -1 / 3 ~ (abc) -1 / 3 , and the 
value 

e = E th (abc) 1/3 = 3 K 

remains constant in time. A Lagrange function of the ellipsoid is written as 

L = Ukin — Up 0 t , Upot = U g + E t h , E t h = i a fo c ^i/3 ' (3.13) 

Equations of motion describing behavior of 3 semiaxes (a, b , c) is obtained from the 
Lagrange function (13.131) in the form 


2>Gm 


du 


5 1 


a = — - 


2 J (a 2 + u) A 3to a ( abc) 3 / 3 
0 


(3.14) 


6 = - 


3 Gm 


du 


+ 


5 1 


2 J (b 2 +u) A 3m b (abc) 1 / 3 
0 


3 Gm 


du 


5 1 


2 J (c 2 + u) A 3 m c (abc) 1 / 3 

0 


(3.15) 

(3.16) 


A 2 = (a 2 + u)(b 2 + u)(c 2 + u). 


4. Dynamic stabilization of non-spherical bodies against unlimited 
collapse, numerical results 

To obtain a numerical solution of equations (13.141) . (13.151) . (13.161) we write them in 
non-dimensional variables. Let us introduce the variables 



t a 


b _ c 

t 

to ’ 

b = 

a 0 ’ a 0 ’ 

m 

P T~r 

U 

Eth - £ 

— . 

, p = —. u = 

— 

Eth — -77—) £ — - 

Uo £0 

ra 0 

Po 

Uo 
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The scaling parameters to, «o, too, po, Uo, £o are connected by the following relations 


t 2 - 

In — 


On „ Gmn Too 


Uo = 


Gmo a o 

System of non-dimensional equations: 


, Po = — 3 -, = Uoao- 


(4.1) 


3m 

a = 2~ a 


du 


+ 


5 1 


(a 2 +u) A 3to a ( abc) 1 / 3 


(4.2) 


b = ~b 


du 


5 1 


0 


[b 2 + u)A 3 m b (abc) 1 / 3 


(4.3) 


3 TO 

c= —c 


du 


+ 


5 1 


0 

,2 — 2 


(c 2 +u) A 3m c (abc ) 1 / 3 


(4.4) 


A 2 = (a 12 + u)(b 2 + u)(c 2 + u). 

In equations (SUD-dUlD only non-dimensional variables are used, and ’’tilde” sign 
is omitted for simplicity in this section. The non-dimensional Hamiltonian (or non- 
dimensional total energy) is: 

H = U kin +U g + E th = ^ (a 2 + b 2 + c 2 )- 


3 m 2 

IF 


du 


1 y ( a 2 + u)(b 2 + u)(c 2 + u) (abc) 1 / 3 


(4.5) 


In case of the sphere (a = b = c, a = b = c) the non-dimensional Hamiltonian and 
non-dimensional equations of motion reduce to: 


H = — ma 2 -to 2 £ 


10 


5a 


(4.6) 


1 


a = — - 


ma* 


TO --£ 


(4.7) 


As follows from dm m for the given mass there is only one equilibrium value of e 

3 to 2 . . 

£eq ^ , ( 4 - 8 ) 

at which the spherical star has zero total energy, and it may have an arbitrary radius. For 
smaller £ < e eq the sphere should contract to singularity, and for £ > e eq there will be a 
total disruption of the star with an expansion to infinity. We solve here numerically the 
equations of motion for a spheroid with a = b 7 / c, which, using gj-flop , are written 
for the oblate spheroid with k = c/a < 1 as 


.. 3 

a = —- 


k- 


arccos k 


2a 2 (l — k 2 )[ y/1 — k 2 \ 3m a (a 2 c) 1 / 3 ’ 


5 1 


c = — 3- 


m 


: 2 (1 - k 2 ) 


1 - 


k arccos k 

VT^k 2 


5 1 


3m c (a 2 c) 1 / 3 ' 


(4.9) 

(4.10) 
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and for the prolate spheroid k = c/a > 1 as 


3 TO 
2 a 2 (k 2 — 



cosh 1 k 
v / F _ l. 


5 1 e 
3 to a (a 2 ^ 1 / 3 ’ 


c = 3 




k cosh 1 k 

y/k 2 — 1 


It is convenient to introduce variables 


5 1 £ 

3to c (a 2 c) 1 / 3 


(4.11) 

(4.12) 


5 £ 

£* = t* = ty/m. 

3 TO Z 


(4.13) 


In these variables the equations (14.9I) - (14.12I) . (14.71) are written as (omitting subscript ”*”) 


a = 


2 a 2 (l — k 2 ) 
3 


c = — - 


k - 


1 - 


arccos k 


y/1 - k 2 \ a (^c) 1 / 3 ’ 


k arccos k 


(i - k 2 ) [ ^/^^F J c (a 2 c y/3 


for the oblate spheroid k = c/a < 1, 

3 


a = — - 


k- 


cosh 1 k 


2 a 2 {k 2 - 1) L yfW^\ \ a (a 2 c)V 3 


1 £ 
H- 


l2 (fc 2 — 1) 


1 - 


A: cosh 1 k 1 

yfW ^I 


c (a 2 c)!/3 


for the prolate spheroid k = cja > 1, and 


(4.14) 

(4.15) 


(4.16) 

(4.17) 


1 — £ 


a = — - 


(4.18) 


for the sphere, where the equilibrium corresponds to s eq = 1. Near the spherical shape 
we should use expansions around k = 1, what leads to equations of motion valid for both 
oblate and prolate cases 


1 — £ 


+ u 


3\ 1 -k 
5 ) a 2 


1 — £ /4£ 4\ 1 — k 

a 2 \ 3 5 ) a 2 

In these variables the total energy is written as 


(4.19) 


and omitting we have 



H, 


H 

9 ’ 

m z 


3 arccos A; 3 £ 

5a y/\ - k 2 + 5 (a 2 c) 1 / 3 ’ 


(oblate) 


(4.20) 



3 cosh 1 k 
5 a y/k 2 - 1 


3 £ 

5 (a 2 c)!/ 3 ’ 


(prolate) 
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1.4- 



time 


Figure 1. Example of regular motion of spheroid with 7 = 4/3, H = —1/5, £ = 2/3. This 
motion corresponds to full line on the Poincare map in Fig|4] 


H = 


10 5 a 

c 2 3 ( 5 2 S 2 

-H-1- 

10 5a V 3 15 


(sphere) 


3 e S 2 S 2 

— H-1- 

5a V 3 9 


(4.21) 


5 = 1 — k, (around the sphere), |<5| <C 1. 

Solution of the system of equations (I4.9I) - (I4. 121) was performed for initial conditions at 
t = 0 : Co = 0 , different values of initial ao, do, fco, and different values of the constant 
parameter e. The results of our numerical calculations are 

(i) sphere (a = b = c): 

e = 1 — total energy equals to zero (H = 0), radius is arbitrary, 
e <1 - spherical star collapses to a singularity, 

£ > 1 - disruption of star with expansion to infinity; 

(ii) spheroid (a = b ^ c): 

at £ ^ 1 the total energy is H ^ 0 - disruption of star with expansion to infinity, 
at e < 1 the total energy is H < 0 - the oscillatory regime is established. In 
the oscillatory regime the dynamical motion prevents the formation of the singularity. 
The type of oscillatory regime depends on initial conditions, and may be represented 
either by regular periodic oscillations, or by chaotic behavior. Examples of two types of 
such oscillations are represented in FigsQ][2] (regular, periodic) and in Fig|3] (chaotic). 
For rigoro us separati on between t hese kinds of oscillations we use a method developed by 
Poi ncare ( Lichtenber g~fc Liebermanl l l983h . The same approach is used in recent papers 
of VandervoortJ ( 2011 . 2014h . 
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Figure 2. Example of regular motion of spheroid with 7 = 4/3, H = —1/5, e = 2/3. This 
motion corresponds to the point inside the regular region on the Poincare map in Fig|4] 
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Figure 3. Example of chaotic motion of spheroid with 7 = 4/3, H = —1/5, e = 2/3. This 
motion corresponds to gray points on the Poincare map in Fig[5] 





































































































































































































































































Dynamics of non-spherical bodies. Zeldovich’s pancakes 


11 



Figure 4. The Poincare map for five regular and two chaotic trajectories in case of 7 = 4/3, 
H = —1/5, e = 2/3. The (a, a) values are taken in the minimum of c. Full black line is the 
bounding curve. The point inside the regular region corresponds to coherent oscillations with 
the same period for a and c values, represented in Fig.2. 


5. Investigation of regular and chaotic oscillations via the Poincare 
sections 


To investigate regular and ch aotic dynamics we use the method of Poincare section 
(j Lichtenber g & Lieber man 19831) , and obtain the Poincare map for different values of the 
total energy H. Let us consider a spheroid with semi-axes a = b 7 ^ c. This system has 
two degrees of freedom. Therefore in this case the phase space is four-dimensional: a, a, 
c, c. If we choose, at given e, a value of the Hamiltonian H 0l we fix a three-dimensional 
energy surface H(a,a,c,c ) = H 0 . During the integration of the equations (ld.9D - (14.12D 
which preserve the constant H , we fix moments ti , when c = 0. At these moments there 
are only two independent values (i.g. a and a), because the value of c is determined 
uniquely from the relation for the hamiltonian at constant H. 

For the same values of H and e we solve equations of motion (14.9I) - (I4.12D at initial 
c = 0, and different a, a. For each integration we put the points on the plane (a, a) 
at the moments ti. These points are the intersection points of the trajectories on the 
three-dimensional energy surface with a two-dimensional plane c = 0, called a Poincare 
section. 


For each fixed combination of e, H we get the Poincare map, represented in Figs0 0 
mm Condition c = 0 splits in two cases, of a minimum and of a maximum of c. The 
Poincare maps are drawing separately, either for the minimum, or for the maximum of c, 
and both maps lead to identical results. The regular oscillations are represented by closed 
lines on the Poincare map, and chaotic behavior fills regions of finite square with dots. 
These regions are separated from the regions of the regular oscillations by separatrix line. 



















12 


G. S. Bisnovatyi-Kogan and O. Yu. Tsupko 




Figure 6. The Poincare map for two chaotic trajectories in case of 7 = 4/3, H = —1/2, 
e = 1/6. The (a, a) values are taken in the minimum of c. Full black line is the bounding curve. 

6. Dynamic stabilization of non-spherical bodies against unlimited 
collapse, discussion 

The main result following from our calculations is the indication to a degenerate nature 
of formation of a singularity in unstable newtonian self-gravitating gaseous bodies. Only 
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Figure 7. The Poincare map for six regular trajectories in the case of 7 = 4/3, H = —3/50, 
e = 9/10. The (a, a) values are taken in the minimum of c. Full black line is the bounding curve. 
The point inside the regular region corresponds to coherent oscillations with the same period 
for a and c values, similar to those represented in Fig[2] 


pure spherical models can collapse to singularity, but any kind of nonsphericity leads to 
nonlinear stabilization of the collapse by a dynamic motion, and formation of regularly or 
chaotically oscillating body. This conclusion is valid for all unstable equations of state, 
namely, for adiabatic with 7 < 4/3. In addition to the case with 7 = 4/3, we have 
calculated the dynamics of the model with 7 = 6/5, and have obtained similar results. 

Note, that region of chaotic behavior on the Poincare map is gradually increasing for 
7 = 4/3 with decreasing of the entropy e and the total energy H. At e = 9/10 and 
H = —3/50 we have found only regular oscillations (FigQ, at e = 2/3 and H = —1/5 
both kind of oscillations are present (Figs|4] [5j) , a nd only chaotic behavior is found 
at e — 1/6 and H = —1/2 (FigEJ. We connect ( Bisnovatvi-Koeran fc Tsunko l2005h 
this chaotic behavior with developm ent of anisotropic instability, when radial velo cities 
strongly exceed the transversal ones ( Antono^ll973t iFridman fo Polvachenko 1985). 

In reality a presence of dissipation leads to damping of these oscillations, and to fi¬ 
nal collapse of nonrotating model, when total energy of the body is negative. In the 
case of core-collapse supernova the main dissipation is due to emission of neutrino. 
The time of the neutrino losses is much larger than the characteristic time of the col¬ 
lapse, so we may expect that the collapse leads to formation of a neutron star where 
nonspherical modes are excited and exist during several seconds after the collapse. In 
addition to the damping due to neutrino emission the shock waves will be generated, 
determining highly variable energy losses during the oscillations. Besides viscosity and 
radiation which can damp the ellipsoid-like motions and allow collapse, there is the pos¬ 
sibility that inertial interactions with higher order modes that must be present in the 
real stratified bodies may cause an inertial cascade and drain the energy from the sec- 
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ond order modes an a non-secular way that does not depend on dissipative coefficients. 
It is very seductive to connect chaotic oscillations with highly variable emissio n ob¬ 


served in the prompt gamma ray emission of cosmic gamma ray bursts (IMazets et al 
198lh . A presence of rotation and ma gnetic fie l d stron gly compli cate the picture o f 


t he core-col l apse s u pernova e xplosion ( Ardelian. Bisnovatvi-Kogan fe Moiseenkol 2005 
1 Moiseenko. Bisnovatvi-Koean fc Ardelian 2005 ). 

In this paper we consider in details only spheroidal bodies. In reality the spheroid will 
become a triaxial ellipsoid during the motion. In addition to spheroids we calculate many 
variants with triaxial figures (see iBisnovatvi-Kogan fc Tsunkol ( 20051) 1. Qualitatively we 
obtain the same results for ellipsoids: no singularity was reached for any e > 0 , and 
establishing oscillatory (regular or chaotic) regime under negative total energy prevents 
the collapse. However, in the case of the ellipsoid with semi-axes a ^ b ^ c we have 
a system with three degrees of freedom and six-dinrensional phase space. Therefore we 
could not carry out rigorous investigation of regular and chaotic types of motion by the 
constructing Poincare map as it was done for spheroid with two degrees of freedom, and 
restricted ourselfs by a description of the spheroidal case. 

In the frame of a general relativity dynamic stabilization against collapse by nonlinear 
nonspherical oscillations cannot be universal. When the size of the body approaches 
gravitational radius no stabilization is possible at any 7 . Nevertheless, the nonlinear 
stabilization may happen at larger radii, so after damping of the oscillations the star 
would collapse to the black hole. Due to development of nonspherical oscillations there is 
a possibility for emission of gravitational waves during the collapse of nonrotating bodies 
with the intensity similar to rotating bodies, or even larger. 


7. Zeldovich’s pancakes and large scale structure of Universe 

One of the most important problem in modern cosmology is a formation of large scale 
structure of the Universe. Observations of larg e scale str ucture of Universe demonstrate 
strong inhomogeneity in density distribution. IZeldovichl ( 1970al lbh has investigated the 
fragmentation of a homogeneous medium under the action of gravitation. It has been 
shown that an arbitrary perturbation in medium with zero pressure leads to formation 
of strongly non-spherical structures (discs), and only pure spherical perturbations lead 
to collapse into the point. These discs are known now as ’Zeldovich’s pancakes’. 

Depending on the medium properties, two types of pancakes can be formed: 

(i) During contraction of the barionic gas the collapse is stopped by the pressure, the 
shock wave is formed. Behind the shock front the cooling, fragmentation and formation 
of vortexes occur, what is necessary for formation of spiral galaxies with large angular 
momentum. 

(ii) During collapse of a non-collisional dark matter the pancake is also formed, but 
no shocks are forming at the stage of maximal contraction. Medium is non-collisional, 
so after initial the collapse, subsequent oscillations occur. The particles fly through each 
other many times, and mixing in the phase space occurs - so call ed v iol ent r elaxation. 
The system reaches quasi-equilibrium state via violent relaxation ( Lvnden-Belll[l967h and 
phase mixing during the most rapid dynamic stages . Established equilibrium state has a 
form, depending on the angular momentum of rotation, and sizes depending on the final 
values of energy and entropy. 

Subsequent nonlinear evolution and interaction between pancakes lead to formation of 
a very complicated large scale structure of matter in the Universe o bserved i n th e sky. 
I t was shown by numerical simulations in two-dimensional space by Doroshkecich_etjy 
( 1980l h The first three-dimensional simulation were performed bv iKlvnin fc Slrandarinl 
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( 1983[). dem onstrating the or igin o f filaments in the dark matter density field. For review 


see 


Shandarin fc Zeldovich ( 1989b where the developments of the original Zeldovich’s 


idea were summarized. The development of perturbations at non-linear stage and struc¬ 
ture formation in the medium with barionic gas and non-collisional cold dark mat- 
ter are investigated by u si ng numeric a l mod e lling by ma n y gro u ps (see, for exam ple. 
Doroshkevich et al ( 19991) . iBoilv et al ( 1999h . Sneth at al ( 2001 1. Boilv et al ( 20021) ). 
For numerical modelling N-body simulations are used, which are very time consuming. 

The first stage of large scale structure formation is a collapse of non-collisional dark 
matter with formation of a pancake and subsequent relaxation with formation of quasi¬ 
equilibrium structures. This stage can be inves tigated using approxima te m odel of large 


scale s tructure formation considered in works of lBisnovatvi-Koganl (|2004r) and lBisnovatvi-Kogan fe Tsupko 
( 20051) . Form of t he collapsing matter hav e been app roximated as rotating homogeneous 
ellipsoid with two ( Bisnovatvi-Koganl2004l) and three ( Bisnovatvi-Kogan fe Tsupkol2005l) 
degrees of freedom. Correct description of pressure effects, attained by such approach, 
and addition of a relaxation permit to get the dynamics of motion without any numerical 
singularities. In our model a motion along three axes takes place in the gravitational Held 
of a uniform density ellipsoid, with account of the isotropic pressure, represented in an 
approximate non-differential way. The relaxation leads to a transformation of the kinetic 
energy of the ordered motion into the kinetic energy of the chaotic motion, and to in¬ 
crease of the effective pressure and thermal energy. All losses are connected with runaway 
particles. The collapse of the rotating three-axis ellipsoid is approximated by a system of 
ordinary differential equations, where the relaxation, and the l osses of energy, mass an d 
angular momentum are taken into account phenomenologically ( Bisnovatvi-Koganll2004j) . 

The system is solved numerically for several parameters, characterizing the configuration. 


8. Equations of motion with dissipation and numerical results 

Equations of motion of rotating dark matter ellipsoid without any dissipation are 
derived in subsection 13.11 see eq. (13.101) . (13.111) . (13.121) . 

In reality there is relaxation in the collisio nless system , connected with phase mix¬ 
ing, which is called ’’violent relaxation” (L vnden-Belll 119671 ). This relaxation leads to a 
transformation of the kinetic energy of the ordered motion into the kinetic energy of the 
chaotic motion and increases effective pressure and thermal energy. 

The main transport process is an effective bulk viscosity. Therefore there is a drag 
force, which is described phenomenologically by adding of the terms 


a 

T~rel 


b 

T~rel 


C 

T~rel 


( 8 . 1 ) 


in the right-hand parts of equations of motion. 

Dissipation leads to a heat production. The rate of a he at production leading t o the 
growth of the entropy of the spheroid matter was found in Bisnovatvi-Kogan i ( 2004t ) from 
equations of motion (13.101) . (13.111) . (13.121) with account of (13.41) . allowing the energy, mass 
and angular momentum losses. 

Here we have scaled the relaxation time T re i by the Jeans characteristic time 


_ 2ti- _ 


27T 


ujj y/AnGp 


= 2-7T 


abc 

3Gm 


( 8 . 2 ) 
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with a constant value of a re i 


T~rel — 2?r OLrel 


abc 

3Gm 


The energy balance relation follows from the dynamic equations in the form 
dUtot dUq dm dU ro t dm i dU ro t dM _ dEfh d£ ^ Ukin Ukin dm 


dt 


dU g dm 
dm dt 


dm dt dM dt 


d£ dt 


T~rel 


m dt 


(8.3) 


(8.4) 


Utot — Ukin T Ug -f- U ro t T Eth. 

The process of relaxation is accompanied also by the energy, mass, and angular mo¬ 
mentum losses from the system. We suggest, that these losses take place only in non¬ 
stationary phases, so these rates are proportional to the kinetic energy Ukin, the charac¬ 
teristic times for mass, angular momentum and energy losses T m i, TMi,T e i are considerably 
greater then T re i. The equations describing different losses may be phenomenologically 
written as 

dUtot U kin dm mUkin TT M / m \ 777 

=- ’ = 77-» U 9< °. TT~ = ( - > (8-5) 

dt 'Tmi dt UgTml M-in \TTlin J 

where mi„ and Mi n are the initial values of corresponding parameters. Scaling all the 
characteristic times by the Jeans value, we have T e i = a e \Tj, T m i = a m iTj, tmi = 
oi-mitj, with constant values of on (i = el, ml, Adi). The function £ determines the en¬ 
tropy of the matter in the spheroid, so that, using m , we have 


dE th | D dU _ 1 d£ ^ _ 4tt u d _ 2 E th 

dt dt ~ (abc) 2 / 3 dt ’ th ~ (abc ) 2 / 3 ’ “ 3 ’ _ 3 V ' 


( 8 . 6 ) 


Using (18.51) in (18.41) we obtain the equation for the entropy function £ in presence of 
different losses in the form 


d£_ 

dt 


(abc) 2 / 3 U k i n 


\( 2 

l 

2 ) 

Urot 

( 2 

1 5 

Ukin 

_ \ 7~rel 

Tel 

1~ml ) 

u g 

\tmi 

T~ml ) 

UgT m i _ 


(8.7) 


The equations, describing approximately the dynamics of the formation of a stationary 
dark matter object include equations of motion (13.101) . (|3.11l) . (13.121) after adding terms 
with relaxation; energy equation (18.71) with account of (13.51) . (13.61) . (13.81) : and equations 
(1531) . describing the losses of the energy, mass and angular momentum. 

To obtain a numerical solution of equations we write them in non-dimensional variables. 
Let us introduce the variables 


t 




b 


b 

— , c 

a 0 


c _ m ~ M _ p 

— , m = — , M = — , p = — , 
a 0 mo M 0 po 


U = ^,E th = 
Uo 

The scaling parameters to, do, mo, Mo, 
relations 


E t h ~ £ Tj 

TT ’ £ j T i ~ , • 

Uo £o t 0 

po, Uo, Uq, £q are connected by the following 


t 2 - 

ir\ — 


Gmo 


M 2 = Gdorng, Uo = 


Gmo 

do 


P o = 


m o 
„ 3 > 


n Mo 2 

Si0 = - O) £0 = tVodn. 

toq ag 


Uo is used for scaling of all kind of energies. Hereinafter the ’tilde’ sign is omitted for 
simplicity. In non-dimensional variables we have m = ^ pabc , M = ^ fl(a 2 + b 2 ), Ti = 

2n an 
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Taking into account violent relaxation, total energy, mass and angular momentum 
losses, the dynamics of the system is described by the following non-dimensional system 
of equations 


a 


OO 

a dm 3to f du 10 1 e 25 M 2 a 

m dt 2 a J {a 2 +u )A + 3 m a {abc) 2 / 3 m 2 {a 2 + b 2 ) 2 T re i ’ 

o 


- —, ( 8 - 8 ) 


■■ b dm 3 m 

m dt 


du 


+ 


10 1 


25 M 2 


2 J (b 2 + u )A 3to 5 {abc) 2 / 3 m 2 {a 2 + b 2 Y T re i 


~ —, (8.9) 


o 


c dm 3 to 
m dt 2 


du 


e = (abc) 2 / 3 U k i 


10 1 


(c 2 + u )A 3m c {abc) 2 / 3 T re i 


1 

2 ) 

U r ot 

( 2 

1 > 

Ukin 

Tel 

7~ml J 

u g 

\tmi 

7~ml ) 

Ug Tmi _ 


( 8 . 10 ) 

(8.H) 


3r m ; 


(a 2 +6 2 +c 2 )/ - 



M 

M~. 


mi, 


iL 

T Ml 


= - , A 2 = {a 2 +u){b 2 +u){c 2 +u). 


( 8 - 12 ) 

This system is solved numerically for s evera l initial parameters . The case of spheroid 
{a = b Y c ) is c onsidered in the paper oflBisnovatvi-Koganl (1200411 . the case of ellipsoid - 
in the paper of Bisnovatyi-Kogan & Tsunkol ( 2005 ). For investigation of the dynamical 
behavior we start the simulation from spherical body of unit mass {mt n = 1), zero or small 
entropy e <C 1. We also specify the parameters, characterizing different dissipations. In all 
calculations with relaxation we use the following values: a re i = 3, a e i = a m i = olmi = 15. 
We set non-zero initial velocities. Because of the rotation around c-axis, the initial sphere 
transforms to spheroid during the collapse, if initial parameters for axes a and b are 
exactly equal {ai n = bi n , aj r = bin). Note that during the motion spheroids may be not 
only oblate but prolate too ( Bisnovatvi-Koganll2004h . To study an appearance of 3-axis 
figures the initial parameters for axes a and b had been slightly different in all variants 
of calculations. 

In the first variant (FigJS]) there is a large initial angular momentum Mi n = 0.5. The 
field of velocities is slightly perturbed by increasing b in comparison with a. Initially we 
observe the collapse and the formation of the pancake. During the motion the difference 
{a—b) increases due to development of a secular instability, and we obtain the transforma¬ 
tion of the Maclaurin spheroid into the Jacobi ellipsoid in the dynamics. Because of the 
relaxation, the oscillating motion damps, and the configuration reaches the equilibrium 
state of the 3-axis ellipsoid. The corresponding behaviors of total energy U to t, entropy 
function e and mass m are represented in Figs [9] [TO] HD We see that the main changes of 
these parameters take place at the stages of first few oscillations: the entropy function e 
increases, the mass m and the total energy Utot decreases. Then these parameters reaches 
the equilibrium values. 

In the second variant of calculations (Fig ll2l) we set a small initial angular momen¬ 
tum Mi n = 0.1. After the collapse, there are temporary oblate and prolate spheroid 
appearance during relaxation. There is no secular instability, and the system reaches the 
equilibrium oblate spheroid of Maclaurin. Here we have an additional instability, there 
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Development of an instability at large angular momentum, and formation of a 
stationary triaxial figure. 



Figure 9. The total energy evolution for the case considered on Fig(8] 
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Figure 10. The entropy function evolution for the case considered on Fig[8] 



Figure 11. The mass evolution for the case considered on Fig[8] 
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Figure 12. The relaxation to the spheroid at small angular momentum. 



Figure 13. Radial instability development in the non-rotating body, without dissipation. 
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Figure 14. Development of the secular instability in a rapidly rotating body, without 

dissipation. 


is increasing of the difference (a — b ) at early stages of motion. However this difference 
remains too small for visualizing in Fig ure. This is the instability characteristic to th e 
system with purely radial trajectories i Antonovl Il973t iFridman fe Polvachenkolll985 1. 
This instability takes place in all variants of calculations, but in most cases it is very 
small and always disappears during the relaxation and formation of the equilibrium fig¬ 
ure. In the cases of large angular momentum this instability does not reveal because the 
difference (a — b) quickly increases as a result of the secular instability. 

For the illustration of the radial instability we consider initial sphere without rotation 
(Mi n = 0) and without any relaxation process. Only three equations (18.811 (18.101) with 
T re i = oo have been used, at constant m = 1, M = 0 and e. We have set small initial 
differences between axes a, b and c, and have obtained an unstable behavior evidently 
visible in FigfTdl 

For numerical investigation of the secular instability we consider equations of motion 
without dissipation, and take the equilibrium spheroids with small perturbation as the 
initial configuration. It is easy to obtain equilibrium parameters of the spheroid from 
equations of motion with zero accelerations. The initial equilibrium configuration of the 
Maclaurin spheroid was slightly perturbed by increasing b. At small angular momentum 
we obtain the conservation of the initial configuration. There are small oscillations around 
the equilibrium spheroid, during which the difference (a — b ) changes sign periodically. At 
large angular momentum we obtain the development of the secular instability, and, after 
several low-amplitude oscillations, the Maclaurin spheroid, is transformed into the Jacobi 
ellipsoid (see Fig lldl where the variant with large angular momentum is represented). 
Because of absence of relaxation, the system does not reach equilibrium configuration and 
oscillations continue. If we include relaxation in this variant, we should observe the final 
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equilibrium state of ellipsoid. Our calculations allow to find approximately numerically 
the point of the onset of the secular instability of the compressible Maclaurin spheroid. 


9. Equilibrium configurations and stability 

Equilibrium of the uniformly rotating figure (spheroid or ellipsoid) is found from the 
equations (18.811 (18.101) with zero time derivatives: 


3 m 

0 = -—a 


0 




du 

+ 

10 1 

6 

25M 2 a 

(9.1) 

(a 2 + u) A 

3to a ( abc ) 2 / 3 

to 2 ( a 2 + b 2 ) 2 ’ 

du 

+ 

10 1 

£ 

25M 2 5 

(9.2) 

(b 2 + u) A 

3to b (abc) 2 / 3 

to 2 (a 2 + 6 2 ) 2 ’ 

OO 

3?7l f 


du 

10 

1 6 

(9.3) 

2 C J 

(c 2 

+ u) A ' 3to 

c (abc) 2 / 3 


Introducing variables 

and integrals 

h{k,ki) = f 
Jo 


u c , b M e 

x=~n, k = ki = j = —, e = —, 
a z a a m m 


(9.4) 


dx 


o (x + 1 ) 3 / 2 \/(fci + x)(k 2 + x) 


I 2 (k, ki) = f -- -— dx — , 

Jo (x + k\) 3 / 2 ^J( 1 + x)(k 2 + x) 

(9.5) 


/ 3 (fc, fci) = f - dx — , 

Jo (x + k 2 ) 3 / 2 ^/(l + x)(kf +a;) 

we write the equations m-m, using (18.121) in the form 

3 to , . 10 e 25j 2 

0--— al 1 (k,k 1 ) + y (^273 + ( i + fc 2 )2 . 

0 _-_ ak 1 I 2 (k,k 1 )+ - ki{kki)V3 + {1 + k 2f 
Excluding to, we obtain 


(9.6) 

(9.7) 

(9.8) 


1 Ii(fc,fci) , x } 15j 2 (fcfci) 2/3 


A : 2 / 3 (fc, fci) 


2e (1 + fc 2 ) 


2^2 ’ 


0 = - 


klh(k,ki) | 1 | 2 15j 2 (fcfci) 2/3 


k 2 I 3 {k,k i) 


+ 1+^1 _ 


26 (1 + fc ?) 5 


. (9.9) 


At a = b, h = 1 equations (19.91) are ide ntical, and determine the equilibrium of Maclaurin 
spheroid, see iBisnovatvi-Koganl ( 2004h . For Jacobi ellipsoids with a ^ b ^ c, we obtain 
the following relation between k and k\ 


J_ 1 J 2 (fc,fci) _ 1 Ji(fc,fci) 

fc 2 A’ 2 / 3 (fc, fci) k 2 h(k,ki) 


= 0. 


E(fc, fci) = 1 


(9.10) 
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This equation has a trivial solution k\ = 1 at all k. corresponding to the Maclairin 
spheroid. Let us find the bifurcation point of the equation (19.101) . at which nontrivial 
solutions appear. While k\ = 1 is always a root of the equation (19.101) . we may write 
F(k, k\) = {k\ — 1 )f(k, fci). Additional root of equation F(k, k\) = 0 appears, when the 
root of the equation f(k,k\ ) = 0 appears at k\ = 1. The root of the zero derivative 
equation F^{k,k\) = f(k,ki) + (ki — l)f' ki (k,ki) = 0 at k\ = 1 coincides with the 
root of the equation f(k, 1 ) = 0 [j], therefore the value of k at the bifurcation point is 
determined by the equation 


dF(k,k i) 

dh 


,=i=0. 


Using m and (19.101) . this equation is written in the form 


(9.11) 


2fc 2 / 3 (fc,fci)| fel= i 


Wi) 

dki 


fci=i 


dh(k,k 1 ) l 
dk\ 


fei=i 


= 0 


(9.12) 


At k\ = 1 we have the analytic expressions 


dh(k,k 0, _ r dh{k, ki) OT r 

|fci=i — --'Oi - ^7 -|fei=i — —oio, 1 o — 


dx 


dk\ 


dk\ 


>o (1 + x) 3 y/k 2 + x 


so that the equation (19.121) is reduced to 

k 2 J 3 = I 0 , where I 3 = I 3 (k, 1 ), and J 2 = h{k, 1 ) = h = h{k, 1 ). 


(9.13) 


(9.14) 


Taking analytically the integrals Iq , It = I 2 and J 3 we have 


(1 - fc 2 )I 0 = -I 2 - - , h = 

1 

h=h = — 


k {1 - fc 2 ) 


1 - 


k arccos k 


VT - k 2 


k(l-k 2 ) 

Using (19.151) in (19.141) we obtain the equation 


k 2 - 


k arccos k 

~7T^W 


(9.15) 


arccosfc fc(13— 10fc 2 ) . a 

v/wF = 3 + 8 fc 2 - 8 fc 4 ’ ) 

which solution k = 0.582724 determines the bifurcation point at the sequence of the 
Maclaurin spheroids. For the uniform spheroid the position of this point does not depend 
on the adiabatic index of the matter. 

Above we have obtained the bifurcation point on the equilibrium curve of the Maclaurin 
spheroids using only the equilibrium relations for the Jacobi ellipsoids. The usual way 
for investigation of stability is to solve linearized equations of motion and to find the 
eigenfrequencies. 

Non-dimensional equations of motion without relaxation are 


OO 
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b = 



du 

10 1 e 

25 M 2 
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( a 2 + u) A 
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f We are grateful to A.I.Neishtadt for useful discussion of this point 






































24 


G. S. Bisnovatyi-Kogan and O. Yu. Tsupko 


3 to f du 

2 C J (c 2 + u)A ' 3 m c ( abc) 2 / 3 
o 


+ , A 2 = (a 2 + u)(b 2 + u)(c 2 + u). (9.19) 


At given a, b, c the equilibrium values of e M (the equilibrium configuration of ellipsoid) 
are obtained using formulae: 


OO 

£ = ^ ° 2 ( a6c ) 2/3 / 


du 


(c 2 + u)A 


(9.20) 


25 M 2 


m 2 (a 2 + b 2 y 


3 to 9 

= —— a 


du 


3 TO 


du 


(a 2 + u)A 


(c 2 + u)A 


(9.21) 


o o 

We linearize equations (19.171) and (19.181) around the equilibrium configuration of the 
Maclaurin spheroid (a = 6, k = c/a). Subtracting (19.171) and (19.181) with using equilib¬ 
rium values (19.201) and (19.211) . we obtain 


6(a — b) = 3 77 
a’ 3 


—Ar 


du 


J (1 + u)(fc 2 + u) 3 / 2 J (1 + u) 3 (k 2 + it) 1 / 2 
0 0 


OO 


du 


8{a-b). (9.22) 


The e, Utot, m and M losses are quadratic to perturbations, so these values remain 
constant in linear approximation. 

Taking 6(a — b) ~ exp(— iut) we come to the characteristic equation 


w 2 = 3 


m 


du 


du 


(1 + u)(fc 2 + u) 3 / 2 J (1 + u) 3 (fc 2 + u) 1 / 2 


(9.23) 


These integrals are expressed in analytical functions: 


w = 


3 m 


k(1 - fc 2 ) 


2arccosfc \ 
(1 — A; 2 ) 3 / 2 / 


3A: 


3 arccos k 


4(1 — A: 2 ) 5 / 2 4(1 — k 2 ) 2 


2(1-k 2 ) 

(9.24) 

For the border of stability w 2 = 0, we obtain the equation again the equation (19.161) . 

The spheroid loses its stability for the transformation into three-axial ellipsoid, at 
the bifurcation point k = 0.58272, e = y/l — c 2 /a 2 = 0.81267. Our approximate equi¬ 
librium equations, even without dissipation, describe the uniformly rotating ellipsoids 
(spheroids), which are not connected by adiabatic relations, and contain a ’’hidden” non¬ 
conservation of the local angular momentum, which preserves the uniformity. Therefore, 
in presence of this ’’hidden” non-conservation, the loss of stability takes place exactly 
in the bifurcation point. The loss of stability is not connected with the relaxation and 
follows from our equations even in the case when all global integrals remain constant. In 
the exact approach the insta bility in this point happens only at the direct presence of 
dissipative terms ( Chandrasekhar! 1969ll . The pure adiabatic spheroidal system preserves 
stability until e = 0.952887, wher e it becomes unstab le via a vibrational mode. 

According to the hypothesis of Ostriker fc Peeble^ (1 9731 ), the stability of an isolated 
axially symmetric system is determined by the ratio U ro t/\U g \. They determined from nu- 
m erical experiments the c ritica l valu e for various co nfiguration as 0.14±0.03. It was found 
in lLai. Rasio fc Shapirol ( 1993 ) and Shapirol ( 20041) that compressible spheroids become 
secularly unstable to triaxial deformations at the bifurcation point, where U ro t/\U g \ = 
0.1375, independent of the adiabatic index n, 7 = 1 + i. Our formula gives exactly the 
same result, which is also confirmed by our numerical simulations. 
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10. Emission of very long gravitational waves and gravitational 
lensing by gravitational waves 

Gravitational radiation is produced during the collapse of the non-sp herical body. 
Gravi tational radiation durin g a formation of a pancake was estimated bv lThuan and Ostrikei 
(1974), and was improved bv INovikovI ( 19751) . who took into account the most important 
stage of the radiation during a bounce. The formula for the estimation of the total energy 
emitted during the collapse and bounce is 


Ugw ~ 0-01 — 


7/2 


*eq 


Me 2 . 


Here a eq is the equilibrium equatorial radius of the pancake, and c„ 


( 10 . 1 ) 


is its minimal 


half-thickness during the bounce, r g = 2 GM/c 2 is the Schwarzschild gravitational radius 
of the body. From our calculations we have a eq /c m i n ^ 100. The value of r g /a eq we 
estimate using the observed average velocity of a galaxy in the rich cluster v g 
km/s, and taking 


3000 


*eq 




10 


-4 


( 10 . 2 ) 


Than the energy carried away by the gravitational wave (GW) may be estimated as 
Ugw ~ 10” 14 M c 2 . If all dark matter had passed through the stage of a pancake forma¬ 
tion, than very long GW with a wavelength of the order of the size of the galactic cluster 
have an average energy density in the universe 


saw ~ 10 PdmC ~ 3 • 10 “ 3 erg/crn . 


(10.3) 


Here we have used for estimation the average dark matter density pd m = 3-10 30 g/cm 3 . 
The average strength Egw of the very long GW may be estimated, taking the relation 

( Landau Lifshidll99,4 

£GM/ = l^G^’ (1 °' 4) 

where h is metric tensor perturbation (non-dimensional), connected with GW, we con¬ 
sider only the scalar, having in mind the averaged value of this perturbation. Taking into 
account h = u>h = 2'Kch/X , A ~ 10 Mpc is the wavelength of GW of the order of the 
size of the cluster. From the comparison of (IIP.3D . (IIP. 41) we obta in the averaged metric 
perturbation due to very long GW in the form ( Bisnovatvi-Koean 20041) 




1/2 


C 2 V 


6-10 


-li 


(10.5) 


for the values of A and £gw-, mentioned above. Incidently we may expect 10 times larger 
amplitude of the GW, than the averaged over the volume value. 

According to general relativity, any gravitational field can change trajectory of pho¬ 
tons or, in other words, deflect light rays. Hence the gravitational field may act as a 
gravitational lens. 

Gravitational lensi ng by gravitational waves in different cases was considered by man; 


uravitational lensing Dy gravitati 
authors (see iBraginskv et al ( 1990i h 


1 by 

Farao nii (1992). lDamour and Esposito-Faresel (199 


and references therein). It was found that the deflection angle vanishes for any lo calized _ 

gravi tational wave packet because of transversality of gravitational waves (Daniour and Esposito-Farese 
1998). Thus if the photon passes throu gh the finite gravitational wave pulse it s deflection 
due to this wave is equal to zero. In Bisnovatvi-Kogan and Tsuokol (20083) we notice 
that the displacement between trajectories of the photon before and after passing the 
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9 Observer 

Figure 15. The observational effect of the displacement in trajectory of the photon. 

wave may occur. On the basis of this result we obtain an approximate formula for the 
estimation of observational effects. 

We confirm analytically vanishing of deflection angle for plane wave pulses. However, 
we have found that the gravitational wave (GW) changes the photon propagation in 
another way, simply shifting its whole trajectory after passing through the GW (see fig. 
USD- This displacement is found analytically for the photon passing through the plane 
GW. Displacement takes place mainly in the case of isolated wave pulses, which have a 
form similar to the top part of the sinusoid or when it has the non-symmetrical top and 
bottom parts of wave profile, and may, in principal, vanish for the periodic wave of a 
long duration. 

Directions of photons passing through the gravitational wave packet does not change, 
therefore any focusing of rays does not occur in this case. Thus the displacement in 
trajectories does not lead to any magnification effect. But the displacement leads to 
change of the angular position of object for distant observer. The change of the angular 


position due to passing of the light ray through the gravitational wave pulse Aa^ is as 
(see fig. [13 



( 10 . 6 ) 


where h is the amplitude of the GW pulse, S is its thickness and D s is a distance between 


the source and the observer. 

Let us estimate the change of the angular position for the GW pulses produces during 
formation of large scale structure of the Universe in dark matter using the following 
parameters: h = 10 —11 , 8 = Mpc , D s = lOOMpc. Than we obtain 


A ad — 2 • 10 8 arcsec. 


(10.7) 






27 


Dynamics of non-spherical bodies. Zeldovich’s pancakes 


Another possible manifestation of gravitational wa ves is influence of them on the po¬ 
larization of cosmic microwave background (lPolnare vlll985f) . It is known that p rimordial 
GWs lead t o formation of B-mod e polar ization ( Seliak fe Zaldarriagalll997 ). see also 
textbook of Gorbunov fe Rubakovl ( 201ll) . This effect is now hunting by different ob¬ 
servational projects, e.g. see BICEP2 (http://bicepkeck.org/). We suppose that B-mode 
polarization formation can be explained not only by primordial gravitational waves but 
also by gravitational waves appeared in the epoch of the formation of the large scale 
structure in the universe. The main distinguishing feature of gravitational waves formed 
during large scale structure formation is anisotropy in direction to rich clusters of galax¬ 


ies. 


11. Conclusions 

We present here the approximate way for investigation of collapse and dynamic sta¬ 
bility based on the model of compressible homogeneous 3-axis ellipsoids. Investigation of 
dynamics is reduced to numerical solution of system of ordinary differential equations of 
semi-axes time evolution. In the presence of dissipation processes there are also equations 
for mass, angular momentum, energy, and entropy evolution. Derivation of equations of 
motion have been performed by variation of the Lagrange function of the ellipsoid. Dif¬ 
ferent equations of state can be easily considered in this model. Right-hand sides of 
equations have analytical form for spheroid case and can be written in terms of ellipti¬ 
cal integrals for 3-axis ellipsoid case. Equilibrium configurations can be also investigated 
by taking zero accelerations in all equations. We apply this method for investigation of 
dynamics of non-spherical stars and dark matter large scale structure formation. 

Dynamic stability of non-spherical stars was investigated. We have solved the equa¬ 
tions that describe, in a simplified way, the Newtonian dynamics of a self-gravitating 
non-rotating spheroidal body after loss of stability. We have found that contraction to a 
singularity occurs only in a pure spherical collapse, and deviations from spherical sym¬ 
metry stop the contraction through the stabilizing action of non-linear non-spherical 
oscillations. A real collapse occurs after damping of the oscillations because of energy 
losses, shock wave formation or viscosity. 

We have investigated the dynamics of 3-axis dark matter ellipsoid. Equations of motion 
for axes of uniform compressible ellipsoid have been obtained by variation of the Lagrange 
function, in which the violent relaxation and losses of a matter, energy and angular 
momentum have been included phenomenologically. 

The system was solved numerically, until the formation of stationary rotating figures 
in presence of the relaxation. For lower angular momentum M we have a formation 
of a compressed spheroid, while at larger M we follow the development of three-axial 
instability and formation of three-axial ellipsoid. The instability in this approximation 
happens at the bifurcation point of the sequence of Maclaurin spheroids, where Jacobi 
ellipsoidal system starts. 

The bifurcation point coinciding with the point of loss of stability is found analytically 
in the form of a simple formula, by static and dynamic approaches. Numerical and an¬ 
alytical considerations give identical results. Development of instability, connected with 
radial orbits, is obtained for slowly rotating collapsing bodies. 

The weak but very long gravitational waves (GW) emitted mainly on the stages of 
the collapse and pancake formation form a long wave GW background. The existence of 
such a long GW in the space between the source and the observer may be registered as 
an action of the gravitational lens. Effect of lensing can be due to displacement of the 
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whole trajectory if light ray is passing through the GW pulse. The displacement leads to 
change of the angular position of object for distant observer. 

The challenging idea is explanation of B-rnode polarization formation not only by 
primordial gravitational waves but also by gravitational waves appeared in the epoch of 
the formation of the large scale structure in the universe. 
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